How far do they go? The outer structure of dark matter halos 
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ABSTRACT 

We study the density profiles of collapsed galaxy-size dark matter halos with masses 10 11 — 
5 • 10 12 Mq focusing mostly on the halo outer regions from the formal virial radius R v i r up to 
5-7R v ir- We find that isolated halos in this mass range extend well beyond R v h exhibiting all 
properties of virialized objects up to 2-3R v i r : relatively smooth density profiles and no systematic 
infall velocities. The dark matter halos in this mass range do not grow as one naively may 
expect through a steady accretion of satellites, i.e., on average there is no mass infall. This 
is strikingly different from more massive halos, which have large infall velocities outside of the 
virial radius. We provide accurate fit for the density profile of these galaxy-size halos. For 
a wide range (0.01 — 2)R v ; r of radii the halo density profiles are fit with the approximation 
p = p s cxp (— 2n[x 1 / n — 1]) + (p m ), where x = r/r s , (p m ) is the mean matter density of the 
Universe, and the index n is in the range n = 6 — 7.5. These profiles do not show a sudden 
change of behavior beyond the virial radius. For larger radii we combine the statistics of the 
initial fluctuations with the spherical collapse model to obtain predictions for the mean and most 
probable density profiles for halos of several masses. The model give excellent results beyond 2-3 
formal virial radii. 
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More than twenty years of extensive work on 
cosmological N-body numerical simulations have 
provided numerous detailed predictions for the 
structure of dark matter halos in the hierarchi- 
cal clustering scenario. Navarro, Frenk, & White 
(1997, hereafter NFW) preceded by the pioneering 
efforts of Quinn et al. (1986); Frenk et al. (1988); 
Dubinski & Carlberg (1991); Warren et al. (1992); 
Crone et al. (1994), suggested a simple fitting for- 
mula to describe the spherically averaged density 
profile of isolated dark matter halos in virial equi- 
librium. Since then numerous simulations were 
done for many relaxed halos of different masses 



and in different cosmologies. The NFW analytical 
density profile 
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has two parameters: the characteristic density p s 
and the radius r s . Instead of these parameters, one 
can use the virial mass of the halo, M v j r , and the 
concentration, C = R v i r /r s . Here the mass M v ; r 
and the corresponding radius R v ; r arc defined as 
the mass and the radius within which the spher- 
ically averaged overdensity is equal to some spe- 
cific value. For the standard cosmological model 
with the cosmological constant ACDM and param- 
eters 51 = 0.3, f^A = 0.7, and h = 0.7 we have 



M v ; r = 47r(340(p m ))i?^ ir /3, where (p m ) is the av- 
erage matter density in the Universe. For a halo 
with this profile, p oc r _1 as r — > and smoothly 
fall off as p oc r~ 3 at the virial radius. The con- 
centration parameter weakly depends on the virial 
mass with a significant scatter comparable to the 
systematic change in C over three decades in M v ; r 
(Bullock et al. 2001; Eke et al. 2001). 

Later simulations paid most of attention to the 
inner slope of the profiles. Some results favored 
a steeper profile than NFW density cusp with 
p oc r -15 (Fukushige & Makino 1997; Moore et 
al. 1998; Jing & Suto 2000; Ghigna et al. 2000). 
More recent simulations of halos with millions of 
particles within R v i r seem to indicate that there 
is a scatter in the inner slope of the density pro- 
files across a wide range of masses - from dwarfs 
to clusters. The inner slope varies between these 
two shapes: the NFW with an asyntotic slope 
of one and the steeper Moore et al. (1998) with 
slope 1.5 (see Klypin et al. 2001; Reed et al. 
2003; Navarro et al. 2004; Dicmand et al. 2004; 
Wambsganss et al. 2004; Tasitsiomi et al. 2004; 
Fukushige et al. 2004). 

Different approximations for density profiles 
were suggested and tested in the literature. Just 
as some other groups, we find that the 3D Scrsic 
three-parameter approximation gives extremely 
good fits for dark matter halos (Navarro et al. 
2004; Merritt et al. 2005). We slightly modify 
this approximation by adding the average matter 
density of the Universe (p m )- This term can be 
neglected, if one fits the density inside the virial 
radius. Yet, at larger distances, it gives an im- 
portant contribution. The approximation can be 
written as 

p{r) = p s exp {-2n[x 1/n - l]j+(p m ), x = r/r s . 

(2) 

where n is the Scrsic index. 

In addition to all the numerical simulations, a 
significant effort has been made to compare the 
predictions of the ACDM model with the observa- 
tions. This is the case of the most recent set of 
high quality observations of large samples of rota- 
tion curves of galaxies or the strong gravitational 
lensing studies which place an important upper 
limit on the amount of dark matter in galaxies 
and clusters in the inner few to tens of kilopar- 
sec (within r s ) where the need of a cuspy density 



profile is still subject of an exciting debate (e.g., 
Florcs & Primack 1994; Moore 1994; de Blok et 
al. 2003; Swaters et al. 2003; Rhee et al. 2003; 
Keeton et al. 1998; Keeton 2001; Broadhurst et 
al. 2004, and references therein). 

On the theoretical side, however, the origin of 
the shape of the dark matter halo density pro- 
file remains poorly understood. It is generally 
accepted that the dark matter halos are assem- 
bled by hierarchical clustering as the result of halo 
merging and continuous accretion. This merging 
scenario has motivated an interest in the analy- 
sis of the mass accretion history of the halos in 
conjunction with their structural properties (e.g., 
Wechsler et al. 2002). The systematic study of the 
NFW density fits to many simulated halos shows 
that their mass accretion history is closely cor- 
related with the concentration parameter C and, 
therefore, with the mass inside the scale radius 
r s (Wechsler et al. 2002; Zhao et al. 2003; Tasit- 
siomi et al. 2004). These results suggest that the 
formation process of the dark matter halos can 
be generally understood by an early phase of fast 
mass accretion and a late phase of slow accretion 
of mass. In this scenario, the inner dense regions of 
the halos are build up early during the fast phase 
of mass accretion when the halo mass increases 
with time much faster than the expansion rate of 
the Universe. At later epochs, during the phase of 
slow mass accretion, the outer regions of the halo 
are built, while its inner regions stay almost intact 
(see Zhao et al. 2003). 

Despite to all this effort dedicated to the un- 
derstanding of the central dense regions of the 
dark matter halos, very little attention has been 
devoted to the study of their outskirts, i.e. the 
regions beyond the formal virial radius. The 
outer parts of the halos and therefore their den- 
sity profiles exhibit in these regions large fluctu- 
ations which can be understood as the result of 
infalling dark matter (including infalling smaller 
halos or substructure) or due to major mergers. In 
both cases the infalling material has not reached 
the equilibrium with the rest of the halo (see 
Fukushige & Makino 2001). On the contrary, 
a considerable observational effort is being made 
to measure the mass distribution around galaxies 
and clusters at large distances using weak gravi- 
tational lensing (e.g., Smith et al. 2001; Guzik & 
Seljak 2002; Kneib et al. 2003; Hoekstra et al. 



2 



2004; Sheldon et al. 2004, and references therein). 
In these cases the distances go well beyond the 
virial radius ranging from few hundred kpc to sev- 
eral Mpc. Individual field galaxies or clusters pro- 
duce a small distortion of the background galaxies 
that allows us to measure the surface mass den- 
sity profile of the dark matter (e.g., Mellier 1999; 
Bartelmann & Schneider 2001). It is customary in 
the weak lensing analysis to specify a dark matter 
halo density profile to model the projected mass 
profile measured with this technique. The NFW 
analytical formula is often adopted and extrapo- 
lated at large distances, beyond R v ; r with p oc r~ 3 . 
This density model may not be accurate enough. 

The motion of satellite galaxies as a test for 
dark matter distribution at large radii (e.g., Zarit- 
sky k White 1994; Zaritsky et al. 1997; Prada et 
al. 2003; Braincrd 2004; Conroy et al. 2004) is 
another observational method, which requires de- 
tailed theoretical predictions for outer density pro- 
files and infall velocities. 

In fact, we really do not know how far the ha- 
los extend. Indeed, the NFW fitting formula was 
proposed and extensively tested to describe dark 
matter halos within R v ; r . This is why the NFW 
density fits are always done within the virial radius 
or even well below the virial radius (< 0.5R v i r ), 
where the halos are expected to be virialized in 
order to avoid such non-equilibrium fluctuations. 
In this context, it is also surprising that the pre- 
diction of the spherical collapse model for the 
mean profile, which can be reasonably expected 
to give good predictions at sufficiently large dis- 
tances (> 2R v ; r ), have not been used. In fact, the 
predictions of this models have only recently been 
worked out (see Barkana 2004) , but they have not 
been tested against numerical simulations. 

Our goal is to carry out a detailed study of the 
density profiles in and around collapsed objects. 
The paper is organized as follows. In Section 2, we 
give the details of the numerical simulations and 
the fits of the density profile of distinct galaxy-size 
halos. In Section 3, we have studied the shape of 
the density and infall velocity profiles of isolated 
halos up to 2 — 3R v ; r for different masses. We 
show in Section 4 that for large radii the mean 
density profile around dark matter halos is in ex- 
cellent agreement with the predictions we have ob- 
tained via the spherical collapse model, which are 
somewhat different from those found by Barkana 



(2004). Finally, in Section 5 discussions and con- 
clusions are given. 

Throughout this paper the formal virial radius 
is the radius within which the mean matter den- 
sity is equal to 340 times the average mean matter 
density p m of the Universe at z = 0. 

2. Numerical simulations and density fits 

The simulations used in this paper are done 
using the Adaptive Refinement Tree (ART) code 
(ART, Kravtsov et al. 1997). The simulations 
were done for the standard ACDM cosmological 
model with Q = 0.3, Cl\ = 0.7, and h = 0.7. 
We study halos selected from four different simu- 
lations, which parameters are listed in Table 1. 

The simulations cover a wide range of scales 
and have different mass and force resolutions. The 
simulation Box20 has the highest resolution, but 
it has only two galaxy-size halos. We use them as 
examples for the structure of halos simulated with 
very high resolution. In most of the cases we limit 
the analysis to well resolved halos: those should 
have more than w20,000 particles inside virial ra- 
dius. The simulation Boxl20 has the larger vol- 
ume, but its mass resolution allows us to use only 
halos with masses larger than 10 13 h~ 1 M Q . Most 
of the analysis of galaxy-size halos is done us- 
ing simulations Box80S and Box80G. In the case 
of the simulation Box80G the whole 80ft-" 1 Mpc 
volume was resolved with equal-mass particles. 
There were about 180,000 halos in the simulation. 
The simulation Box80S was done using particles 
with different masses. Only a small fraction - a 
10/i -1 Mpc radius region - of the box was resolved 
with small-mass particles. The high resolution re- 
gion has an average density about equal to the 
mean density of the Universe. It was chosen in 
such a way that the halo mass function in the re- 
gion is representative for the typical region of this 
size. For example, the region does not have a mas- 
sive cluster. The most massive halo in the region 
has mass 2.6 x 10 13 /i -1 M©. There are 5 halos with 
mass larger than 10 13 ft. _1 M©. Altogether, there 
are about 60,000 halos in this simulation. 

The density profile of each halo, which we study, 
is fit by the approximation given in eq.( 2). Each 
fit provides the concentration C and the Sersic in- 
dex n. We often average density profiles of halos 
for some range of M v i r . When doing so, we scale 
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the radii to units of the virial radius of each halo 
and then average the densities. The averaged den- 
sity profile is then fit again. In some cases, before 
we do the fitting, we also split the halo population 
of a given mass range into 3-4 sub-samples with 
a narrow range of concentrations. The parame- 
ters of the density fits together with some other 
properties of the halos are given in Table 2. 

The halos in our catalogs come from different 
environments. Some of them are inside the virial 
radii of larger halos; some have strong interac- 
tions with smaller, but still massive neighbors. We 
call a halo "distinct" if it does not belong to a 
larger halo. Most of the time we are interested in 
isolated halos: distinct halos, which do not have 
large companions. We search for halos around the 
given halo. If within the distance d x R v j r the 
largest companion is smaller than M v j r /m, then 
the halo is called isolated. When doing the pair- 
wise comparisons of the halos, we use the largest 
virial radius of the two halos, but we use M v ; r of 
the given halo for the test of the masses. Different 
isolation criteria are used. We typically use the 
d = 2 and m = 5 combination (no massive com- 
panion within 2R v ; r ). For Milky- Way size halos 
with M v ; r w 10 12 /i _1 Mq this condition typically 
gives 50%-60% of all distinct halos of this mass. 

3. Halo profiles and infall velocities 

Figure 1 gives examples of density profiles of 
two halos with virial masses 1.4 x 10 12 /i~ 1 M 
(left panel) and 2.6 x lO 11 /^ 1 !/© (right panel) in 
the simulation Box20. The halos have the virial 
radii of 230/i~ 1 kpc and 130/i _1 kpc respectively. 
The halos were done with very high resolution, 
which allows us to track the density profile be- 
low 0.01R v i r . The larger halo is isolated with its 
nearest companion being at 3.5R v ; r . The density 
profile of the halo has some spikes due to substruc- 
ture, but otherwise it clearly extends up to 3R V h- 
where we see large fluctuations due to its compan- 
ion. The smaller halo on the right panel has a 
neighbor at 2R v ; r . So, it is not isolated. Eq.( 2) 
gives very good approximations for both halos. 

Figure 2 gives more information on the struc- 
ture of these two halos. The 3D rms velocities are 
what one would naively expect for "normal" halos. 
The rms velocity first increases when we go from 
the center and reaches a maximum at some dis- 
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Fig. 1 . — Dark matter density profiles of two dark 
matter halos (full curves) in the simulation Box20. 
The halos have virial masses of 1.4 x 1O 12 /i _1 M0 
(left panel) and 2.6 x 10 11 h,~ 1 M Q (right panel). 
The larger halo has a neighbour at 3.5 R v ; r which 
is the halo on the right panel. This smaller halo is 
responsible for the spike at large radii in the den- 
sity profile. In turn, the halo on the right panel 
has its own smaller neigbour at 2R vir observed as 
a spike and an extended bump in the density pro- 
file. The dashed curves show the 3D Sersic profiles. 
The halo density profiles extend well beyond the 
formal virial radius with the Sersic profile provid- 
ing remarkably good fits. 

tance. The radius of the maximum rms velocity is 
smaller for the halo on the right panel. This is be- 
cause it has larger concentration. At larger radii 
the rms velocity first declines relatively smoothly. 
It has some fluctuations due to substructure. At 
radii larger than R v ; r the decline stops. The aver- 
age radial velocity is more interesting and to some 
degree is surprising. Nothing unusual inside R v i r : 
it is practically zero with tiny (ss 5km s _1 ) varia- 
tions due to substructure. This is a clear sign of a 
virialized object. At larger distances the fluctua- 
tions in the velocity increase, but not dramatically 
if we compare those fluctuations with the rms ve- 
locities. The smaller halo on the right has a nar- 
row dip (V ta d ~ —20km s _1 ) at 1.2R v ; r apparently 
due to a satellite, which is moving into the halo. 
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Table 1 : Parameters of Simulations 



Name 


Box 


Mass resolution 


Force resolution 


Number of particles 


Number of steps 




(ft" 1 Mpc) 


(r'M ) 


(h,- 1 kpc) 


(io 6 ) 




Box20 


20 


6.1 x 10 5 


0.15 


9.0 


5 x 10 5 


Box80S 


80 


4.9 x 10 6 


0.15 


160 


1 x 10 6 


Box80G 


80 


3.2 x 10 8 


1.2 


134 


1.25 x 10 5 


Boxl20 


120 


1.1 x 10 9 


1.8 


134 


1.25 x 10 5 



Table 2: Parameters of halos 



Model Name 


Virial Mass 


Number of halos 


Number of particles 


Concentration 


n 


Figure 




(/•"'Me) 




inside _R vir 








Box20 


1.4 x 10 12 


1 


2.2 x 10 6 


14.3 


5.9 


1 




2.8 x 10 11 


1 


4.6 x 10 5 


16.7 


7.5 


1 


Box80S 


(1 ±0.3) x 10 11 


162 


(2 ± 0.7) x 10 4 


18 


7.0 


3 




(4± 2) x 10 11 


79 


(8 ± 4) x 10 4 


14.0 


6.5 


3,4 


Box80G 


(6 ±2) x 10 12 


192 


(1.8 ±0.6) x 10 4 


11.7 


6.3 


3 



The surprising result is what we do not find, i.e., 
there is no infall on the two halos. This may seem 
like a fluke. Indeed, halos must grow. Their mass 
must increase with time. In order for the mass to 
increase, there should be on average negative in- 
fall velocities just outside the virial radius. Yet, 
we do not find those. We will later see that these 
two halos are not flukes, but are typical examples 
for halos of this mass range. What is important at 
this stage is that the infall velocities outside of R v i r 
are very small. They are significantly smaller than 
the rms velocities. These small radial velocities in- 
dicate that halos may extend to radii significantly 
larger than their formal virial radii. 

The average density profiles of isolated halos of 
different masses are shown in Figure 3. The isola- 
tion criteria in this case is no massive satellite in- 
side 2R v j r (d = 2, m — 5; see Sec 2). The smooth 
density profiles, which are extremely accurately fit 
by eq.( 2), extend from the smallest resolved ra- 
dius all the way to 2R v j r . Parameters of the den- 
sity fits are given in Table 2. Note that in order to 
reduce the range of variations along the y-axis, we 



plot density multiplied by (r/R vir ) 2 . The horizon- 
tal parts of the curves in this plots correspond to 
density declining as p oc r~ 2 . The density profiles 
are well above the average density of the Universe 
throughout all the radii. Even at 5R v i r the aver- 
age density profile is still 4-5 times larger than the 
mean density. The upturn at large radii tells us 
that the density declines less steep than r~ 2 . 

The NFW approximation provides less accurate 
fits as shown in Figure 4. One may attribute the 
success of the 3D Sersic approximation to the fact 
that it has three free parameters, while the NFW 
has only two. This is correct only to some degree. 
The problem with the NFW is that it has a slightly 
wrong shape at radii around r s : its curvature is a 
bit too large. In Figure 4 this is manifested by 
an extended hump close to the maximum of the 
curves at r w (0.05 — 0.2)R v ; r . One can shift the 
NFW slightly to the right and down to make the 
fit more accurate for most of the body of the halo 
(r > 0.05R vir ). In this case the NFW fit goes be- 
low the halo density at small radii r < 0.05R v i r , 
which sometimes was interpreted as if the central 
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Fig. 2. — The top curves on each panel show 
the 3D rms velocities at different radii for the ha- 
los shown in Figure 1. The bottom curves show 
the average radial velocities. The radial velocities 
show that for isolated halos of Milky- Way size and 
smaller halos there is no infall velocities outside of 
formal virial radius. 

slope is steeper than -1. Overall, the 3D Sersic ap- 
proximation provides remarkable accurate fit. For 
r = (0.01 — 2)R v j r the errors are smaller than 5%. 

Figure 5 shows how the density profile depends 
on particular choice of the isolation criterion. In 
this case we selected few hundred halos in the sim- 
ulation Box80G with masses M « 10 12 h~ 1 M Q . 
Qualitatively the same results are found for halos 
with different masses. Conclusions are clear: More 
strict isolation conditions result in smaller density 
in the outer parts of halos with almost no effect 
inside the virial radius. Even outside of the virial 
radius the difference are not that large. The dif- 
ference between distinct (not isolated) halos and 
halos, which have no massive companions inside 
2R vir , are not more than a factor ~ 1.5. To large 
degree, this is not surprising because our isolated 
halos are typical, i.e., more that 1/2 of the halos 
in this mass range are "isolated" . 

One of the misconceptions, which we had be- 
fore starting the analysis of the outer regions of 
dark matter halos is that at large distances the 




Fig. 3. — Average density profiles for halos with 
different virial masses. The 3D Sersic profile pro- 
vides very good fit with few percent errors within 
2R V ir- Even at 3R V ir the error is less than 20-30 
percent. The density profiles are well above the 
average density of the Universe throughout all the 
radii. 



deviations from halo to halo are so large that it 
is very difficult to talk about average profile or a 
profile altogether. This appears to be not true. 
Figure 6 shows the halo to halo rms deviations 
from the average density profile for halos of a given 
mass. In order to construct the plot, we split 
the halo population into three ranges of concen- 
trations and found the average and deviations for 
each concentration bin. Then the results of differ- 
ent concentrations were averaged. This splitting 
into concentrations is needed only for the central 
region r < 0.1R V i r because here the average profile 
depends on the concentration. This plot demon- 
strates that there is no drastic change in the devi- 
ations at the virial radius. The deviations increase 
with the distance, but they are not unreasonable. 

Figure 7 shows the average radial velocity pro- 
file for halos of vastly different masses. We used 
many dozens of halos for each mass range. Just as 
in the case of the two individual halos in Figure 1 
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Fig. 4. — Comparison of the NFW and the 3D 
Sersic fits. The thick full curve in the bottom 
panel shows the average density for halos with 
mass < M >= 3 x 10 11 h- 1 M Q . The 3D Sersic 
fit is the dashed curve. The NFW fit is presented 
by the thin full curve. The top panel shows the 
errors of the fits. The full curve is for the NFW 
approximation, and the dashed curve is for the 
3D Sersic fit. In outer regions r = (0.3 — 2)R vir 
both fits have practically the same accuracy. Both 
fits start to fail at larger distances. Overall, the 
3D Sersic approximation provides remarkable ac- 
curate fit. 

and 2, there is no systematic infall of material be- 
yond formal virial radius for small virial masses. 
The situation is different for group- and cluster- 
sized halos (two top panels). For these large halos 
there are large infall velocities, which amplitude 
increases with halo mass. 

4. Predictions from the spherical collapse 
model 

Here we obtain the predictions of the spheri- 
cal collapse model (Gunn & Gott 1972) for the 
mean halo density profile and compare them with 
the results found in the numerical simulations. 
The spherical collapse model provides a relation- 
ship between the present nonlinear enclosed den- 
sity contrast, 5 = (p(r))/(pm) — F within a sphere 




r/R, 



Fig. 5. — Effects of different isolation criteria on 
the density profiles. The dashed curve shows the 
average density profile of distinct halos with mass 
M w 1O 12 /i _1 M in Box80G. The dot-dashed and 
full curves are for isolated halos, which do not have 
massive companions inside 2R v j r and 3R v i r respec- 
tively. More strict isolation conditions result in 
lower density outside of the virial radius. At all 
shown radii the densities are well above the aver- 
age density of the Universe. 

of given radius r and the enclosed linear density 
contrast, Si, (i.e., the initial fluctuation extrapo- 
lated to the present using the linear theory) within 
the same sphere. This relationship along with the 
statistics of the initial fluctuations furnishes defi- 
nite predictions for the mean density profile. 

We expect the spherical collapse model to give 
the mean density profile around virialized halos for 
sufficiently large radii - significantly larger than 
R v i r . In fact, for r/R v i r > 2.5 we find a good 
agreement between the prediction of the spherical 
model and the mean density profiles obtained from 
numerical simulations for several mass ranges. By 
mean density profile we imply any representative 
profile for a given mass range. We consider three 
different representative profiles: the most proba- 
ble profile, the mean profile and, what we call the 
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Fig. 6. — RMS deviations of halo profiles from av- 
erage profile for halos of given mass. The curves 
in the plot show results for halos with masses indi- 
cated in the plot. The deviations inside the formal 
virial radius are clearly smaller than for the outer 
regions. Inside the virial radius the deviations are 
smaller for halos with smaller mass. 

typical profile. This last profile is simply the mean 
profile in the initial conditions spherically evolved. 

Our procedure to obtain the typical density 
profile will be presented in full detail in Betancort- 
Rijo et al. (2005). In essence, we use the statis- 
tics of the initial Gaussian field to obtain the 
mean profile of the enclosed linear density con- 
trast, Si(q), around a proto-halo as a function 
of the Lagrangian distance to the center, q. We 
then use the spherical collapse model to obtain 
the present density contrast, 5(r), as a function of 
the present comoving (i.e., Eulerian) distance to 
the center of the halo, r. 

If Q is the virial radius in Lagrangian coordi- 
nates, then the enclosed linear contrast 5i(q,Q) 
must satisfy the condition that 5i(q = Q) = S V i r , 
where S V i r is the linear density contrast within the 
virial radius R v i r at the moment of virialization. 
The present average enclosed fractional density 
5(r = R v i r ) is equal to 340 (the value of over- 
density A v ; r used in our simulations to define the 
virial radius). To obtain 5 V i r we must find the <5; 
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Fig. 7. — Average radial velocities for halos with 
different virial masses. The full curves show re- 
sults for isolated halos and the dot-dashed curves 
are for distinct halos. The velocities are practi- 
cally zero within 2 — 3 R v i r for halos with mass 
smaller than 1O 12 /i _1 M . The situation is differ- 
ent for group- and cluster- sized halos (two top 
panels). For these large halos there are large in- 
fall velocities, which amplitude increases with halo 
mass. 

corresponding to a present density contrast equal 
to 340. This condition leads to a S V i r value of 1.9. 
We shall later comment on this once we introduce 
the spherical collapse model. 

This condition simply ensures that the proto- 
halo evolves into an object which at present is 
virialized within the prescribed virial radius R v ir- 
If this were the only constraint on the linear pro- 
file (equal to the initial profile except for an overall 
factor), it have been shown (see Section 4 in Patiri 
et al. 2004) that: 



5i(q, Q) ~ S mr exp [ - b((^) 2 - l)] = 5 {q, Q), 



(3) 
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where b is a coefficient depending on R v i r . ft is 
given below. 

However, although this profile does not lead to 
wrong results, in order to achieve the accuracy 
required here and to obtain the correct dependence 
of the shape on mass of the halos we must include 
an additional constraint. This constraint is: 



8i(q, Q) < 8 vir V q > Q. 



(4) 



It means that there are no radii larger than R v i r , 
where the density contrast is 340 or larger. Oth- 
erwise, the virial radius would be larger than R v i r - 
In the large mass limit this condition becomes ir- 
relevant, and the linear density contrast is given 
by eq.(3). However, in the general case the mean 
linear profile, 6i(q,Q), is somewhat steeper than 
8o{q,Q)- 



a'(q, Q) exp[-x 2 



8i{q,Q) = 8 (q,Q) 

Svir ~ $o(q, Q) 



r;erfc[— x] 



x = 



V2a'(q,Q) 
with So(q,Q) given in cq.(3), and 



(5) 
(6) 



where a(q) and a(Q) are the rms linear density 
fluctuations in spheres with lagrangian radii q and 
Q. We use as = 0.9 as in the simulations. 

With these definitions we have for the previ- 
ously defined function 6(R v j r ): 



6(i?vir) 



1 dlncr(Q) 
I dlnQ 



Q=R vir (340) V3 



(8) 



Using the linear density contrast eq.(5-6) we 
can now obtain the nonlinear density contrast, 
S(r). For any given radius r the nonlinear den- 
sity contrast S(r) is given by the solution to the 
equation: 

6 l (6(r))=6 l (q,Q), (9) 
q = r(l + S(r)) 1 / 3 , Q = i? OTr (340) 1/3 . (10) 



The derivation of this eq.(9) may be found in 
Patiri et al.(2004), although here the equation is 
presented in a slightly different form. The left 
hand side of this equation is simply the relation- 
ship, Si (8), between the present and the linear S 
value in the spherical collapse model (see Shcth 
& Tormen 2002) evaluated at 8 = S(r). For the 
cosmology considered here we have for 8i(8): 



8i(8) 



1.68647 



1.35 



1.676 
1.68647 
1.12431 

(1 + S) 1 / 2 + (1 + 5)0.58661 



(1 + 5) 2 / 3 

0.78785 



(11) 



The right hand side of eq.(9), which depends on 
5(r) through q, is given by expression (3). 

Inserting eq.(5-6) into eq.(9) with 5 v i r = 1.9, 
using for b the values 0.186 and 0.254 for the two 
masses 6.5 x lO^/i" 1 ]^© and 3 x 10 12 h^MQ dis- 
cussed here, and solving for 8(r), we obtain the 
profiles for the present enclosed density contrast, 
S(r). To obtain the density S'(r) at a given radius 
r we simply need to take a derivative: 



( 1+ , f(r)-I^* W . 

(12) 

where (p m ) is the mean matter density of the Uni- 
verse. 

We must now comment on the value of 8 v i r that 
we use. If the standard spherical collapse model 
were valid down to the virial radius we could use 
expression (11) to obtain it: 

S vir = <5z(340) = 1.614 

However, we know that this is not true because 
bellow ~ 2 virial radius there is substantial 
amount of shell-crossing that render the men- 
tioned model unappropiate. This causes 8 to grow 
much more slowly with Si , so that when S takes the 
value 340, Si takes the value 1.9 (Betancort-Rijo et 
al. 2005). 

To obtain the predictions for the most proba- 
ble and mean profiles, the probability distribution, 
P(8, s), for the value of S at a given dimensionless 
radius s = r/R vir is needed. This distribution is 
derived in Betancort-Rijo et al.(2005) along the 
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same lines we have described for the the typical 
profile. Here we simply use it to compute the most 
probable, 5(s) pro b profile which is is simply given 
by the value at which P(S, s) has its maxima. The 
mean profile, < 5(s) > is obtained from: 

/oo 
P(5,s)dd. (13) 

In practice, since the approximation we use for 
P(S, s) has an unduly long tail, we have to trun- 
cate it artificially. So, we only integrate up to 
a S value where P(S, s) has fallen to a twenty 
fifth of its maximum value. The results are pre- 
sented in Figure 8, where we have computed, both 
for S and for 5', the most probable (squares) and 
mean (crosses) density profiles found in our simu- 
lations for two mass intervals with the mean val- 
ues equal to the masses used in the theoretical 
derivation. We have taken 277 halos in the mass 
range (6.5 ± 1.5) x 10 w h- 1 M Q from BoxSOS and 
654 halos in the mass range (3± 1) x 10 12 h~ 1 M Q 
from Box80G. No isolation criteria was used. 
In Table 3 we list for the mean halo with mass 
< M >= 3 x 10 12 /i _1 M Q the estimations of the 
most probable and mean value of the density at 
different radii compare with that from the spheri- 
cal collapse model for the most probable, the mean 
and the typical profiles. 

In Figure 8 one can see that beyond 2 virial 
radius the mean and the most probable profiles, 
both for S and for 5', differ considerably. This is 
due to the fact that for this radii the probability 
distribution for S, P(5,s), is rather wide, with a 
long upper tail. This can be seen in Figure 9 were 
this distribution is shown inside 3.5 ± 0.05 virial 
radius for the mass < M >= 3 x lO 12 /i _1 M . We 
show for comparison the theoretical prediction for 
P(S, s) as well as we give the most probable <5 max 
and mean value < 5 > of the distribution. 

It is apparent from Figure 8 that the 5' pro- 
files are steeper for smaller masses, so that they 
go below the background at smaller r/R vir and 
reach larger underdensities. The S profiles are also 
steeper for smaller masses although the difference 
is, obviously, much smaller. We have found that 
the theoretical prediction for the typical and the 
most probable profile arc, in general, almost indis- 
tinguishable (see Table 3). They are both found 
to be in very good agreement with the most prob- 
able 5 profile found in the numerical simulations 



beyond two virial radii. There is also qualitative 
agreement between the predictions of the most 
probable 5' profiles and those found in the numer- 
ical simulations. It must be noted, however, that 
by predicted 8' profile we understand simply the 
one obtained from the corresponding S profile by 
means of relationship given in cq.(12). Note that 
this is not the same as the most probable profile 
for 5' (see Figure 8). This is due to the fact that 
the most probable 5' value at a given s(= r/R V i r ) 
corresponds to a different halo than the most prob- 
able 5 value at the same s. This explains that, 
while the prediction for 5 pro b agrees very well with 
the simulations, the agreement is not so good for 
dprob- Tne d'prot, obtained from S prob by means of 
expression eq.(12) is not a proper prediction but 
an indicative value, since we can not envisage a 
feasible procedure to obtain a proper prediction. 
On the contrary, the mean profiles < S > are ex- 
actly related to < 6' > by means of expression 
eq.(12). The predictions for both profiles arc in 
good agreement with the numerical simulations, 
showing a much flatter profile beyond 2 virial ra- 
dius than those corresponding to <5 pro b and <5 prob . 
This agreement is remarkable given the fact that 
the expression used for P(S, s) is only a first ap- 
proximation (see Betancort-Rijo ct al. 2005). 

It is interesting to note, as we have previously 
pointed out, that larger masses have somewhat 
shallower profiles. In order to predict this trend 
correctly we must use the initial profile given by 
eq.(5-6). If we dropped the second constraint (that 
is given in eq.(4) and use in eq.(9) the initial profile 
given by eq.(3), which corresponds to high mass 
objects, the prediction would be the opposite. The 
reason for this being that, in this limit, the initial 
profile depend on mass only through c which in- 
creases with increasing mass, thereby leading to 
steeper profiles for larger masses. 

As we stated in the introduction, the computa- 
tions of the mean profiles (in fact, the typical pro- 
file) around halos has independently been made 
by Barkana (2004). He used the spherical model 
and, in principle, imposed on the initial profile the 
same constraint as we do. The computing proce- 
dure he followed was somewhat different involving 
some approximations. He do not give explicitly 
the equation defining the profile, so that accu- 
rate comparison with our results are not possible. 
Furthermore he used different values of S V i r , A vir . 
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Fig. 8. — The most probable (squares) and the mean (crosses) halo density profiles up to 7R v i r for the two 
masses (< M >= 3 x 10 12 ft- _1 M Q and < M >= 6.5 x 10 10 /i. -1 Mg in our simulations. We compare the 
simulated data with the predictions from the spherical collapse model for the most probable (solid line) and 
typical profiles (dashed line) which are almost indistinguishable and the mean profiles (dot-short dashed 
line). We have estimated that the errors in the most probable values of S and 5' are about 25% larger than 
the error showed in the plots for their mean values. 



However, his results are in good qualitative agree- 
ment with our predictions for the typical profile. 

We have so far considered randomly chosen ha- 
los, i.e. without isolation criteria. When halos are 
chosen according with an isolation criteria they 



differ from the randomly chosen ones in two re- 
spects. On the one hand, the probability distri- 
bution for S at a given value of s is narrower, so 
that most profiles cluster around the most prob- 
able one. Therefore, the difference between the 
mean and the most probable density profile be- 



ll 



Table 3: Comparison between the simulated mean halo density profile and the theoretical predictions from 
the spherical collapse model for the mass < M >= 3 x 10 12 h^ 1 M Q . 



Numerical Simulations Spherical Collapse 



r/R v ir 


< <5 > 


^prob 


<<5' > 


si 
°prob 


< <5 > 


^prob 


St 


<S'> 


St 

°prob 


Si 


1.0 


337 


323 


69.1 


52.4 


531 


398 


385 


40.4 


2.5 


7.8 


1.5 


129 


116 


30.8 


15.1 


184 


115 


118 


22.6 


1.9 


4.3 


2.0 


67.7 


54.6 


19.5 


3.73 


84.8 


53.7 


51.7 


14.7 


1.5 


2.8 


2.5 


43.2 


31.1 


15.4 


1.6 


49.3 


28.6 


27.6 


10.3 


1.2 


1.9 


3.0 


30.8 


17.4 


11.2 


0.78 


32.2 


17.3 


16.6 


7.6 


1.0 


1.4 


3.5 


22.1 


12.6 


8.5 


0.58 


22.6 


11.3 


10.9 


5.8 


0.79 


0.97 


4.0 


18.4 


9.0 


5.9 


0.18 


16.8 


7.9 


7.6 


4.5 


0.61 


0.70 


4.5 


14.1 


5.8 


5.5 


0.09 


13.0 


5.8 


5.5 


3.5 


0.46 


0.50 


5.0 


12.2 


4.4 


4.3 


-0.08 


10.4 


4.3 


4.1 


2.7 


0.32 


0.34 


5.5 


9.4 


3.5 


3.3 


-0.16 


8.5 


3.3 


3.2 


2.1 


0.20 


0.22 


6.0 


7.9 


2.8 


3.1 


-0.17 


7.1 


2.5 


2.5 


1.6 


0.09 


0.12 


6.5 


6.8 


2.1 


2.5 


-0.22 


5.7 


2.0 


2.0 


1.2 


-0.008 


0.03 


7.0 


6.2 


1.6 


1.9 


-0.33 


4.6 


1.6 


1.6 


0.9 


-0.10 


-0.04 



Note. — The symbols < S >, <5 pro f,, St stand, respectively, for mean, most probable, and typical averaged enclosed fractional 
ovcrdensity. The corresponding primed simbols are for the local fractional overdensities at given radius. 



comes smaller. On the other hand, isolated pro- 
files lay, on average, on somewhat more under- 
dense environment than non-isolated ones, so that 
their most probable profiles are slightly steeper. 
Both effects may be seen by comparing Table 3 
and Table 4 or the upper right pannel in Figure 8 
with Figure 10 where we show the local density 
profile for the isolated mean halo density profile 
for the mass range 3±1 x lO 12 /i~ 1 Af . In this mass 
range we have selected halos that do not have a 
companion with mass larger than 10% of the halo 
mass within 4R v j r . In total there are 156 halos, 
i.e. one quarter of all halos in this mass range. 

5. Conclusions and Discussions 

We perform a detailed study of the density pro- 
files of isolated galaxy-size dark matter halos in 
high resolution cosmological simulations. We de- 
vote careful consideration mainly to the halo outer 
structure beyond the formal virial radius R v i r . We 
find that the 3D Sersic three-parameter approx- 
imation provides excellent good density fits for 
these dark matter halos up to 2-3R v ; r . These pro- 
files do not display an abrupt change of shape be- 



yond the virial radius. The halo-to-halo rms devi- 
ations from the average profile for halos of a given 
mass show that there is no a drastic change in the 
deviations at the virial radius. We show that these 
density profiles differ considerably from the NFW 
density profile beyond 3R v i r where the density pro- 
file are slower than r~ 2 . This result must not be 
seen as a contradiction when is compared with the 
r -3 NFW fall-off at large radii since we must re- 
member that the NFW analytical formula was pro- 
posed and extensively tested to describe the struc- 
ture of virialized dark matter halos within R v i r - 
Although surprising, it is customary, for example, 
to see that the weak lensing analysis is often done 
with the NFW fit extrapolated to distances well 
beyond R v j r , up to large distances of several virial 
radii (up to few Mpc). This approach may not be 
accurate enough given the results presented in this 
work. 

We also find that the isolated galaxy-size ha- 
los display all the properties of relaxed objects up 
to 2-3R v i r . In addition to their relatively smooth 
density profiles seen at large radii, by studying 
halos average radial velocities, we find that there 
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r/R, 



Fig. 9. — Distribution of the fractional cumulative 
density S inside 3.5 ± 0.05R v ir for the mean halo 
of mass < M >= 3 x 10 12 /i -1 M©. We show for 
comparison the theoretical prediction of P(S, s) as 
well as we give the most probable 5 max and mean 
value < S > of the distribution. We display the 
density distribution from its minimum value up to 
la from its mean < 5 > (8% of the values of the 
distribution are beyond la). 

is no indication of systematic infall of material 
beyond the formal virial radius. The dark mat- 
ter halos in this mass range do not grow as one 
naively may expect through a steady accretion of 
satellites, i.e., on average there is no mass infall. 
This is strikingly different for more massive ha- 
los, such as group- and cluster-sized halos which 
exhibit large infall velocities outside of the formal 
virial radius. For large halos the amplitude of the 
infall velocities increases with halo mass. 

For larger radii beyond 2-3 formal virial radius 
we combine the statistics of the initial fluctuations 
with the spherical collapse model to obtain pre- 
dictions of the mean halo density profiles for ha- 
los with different masses. We consider two pos- 
sibilities: the most probable and the mean den- 
sity profiles. We find that the most probable pro- 
file obtained from our simulations is in excellent 
agreement with the predictions from the spheri- 
cal collapse model beyond 2-3 virial radius. For 



Fig. 10. — The mean (crosses) and most prob- 
able values (squares) of the local density profile 
for the mean isolated dark matter halo with mass 
< M >= 3 x lCP^M©. 

the mean density profile the predictions are not so 
accurate. This is due to the fact that the approxi- 
mation, which we are using for the distribution of 
5 at a given radius, has an artificially long tail (we 
are presently working on a better approximation). 
Even so, the predictions are qualitatively good and 
quantitatively quite acceptable. We think that 
the discrepancies between the data and the pre- 
dictions at radii smaller than 2-3 virial radii are 
due to the fact that these inner shells are affected 
by the shell-crossing. We hope that an appropiate 
treatment of this circunstance will lead to accurate 
predictions for all radii, so that the mean spheri- 
cally averaged profiles may be understood in terms 
of the spherical collapse. We find the results pre- 
sented here very encouraging in this respect and 
we are currently working along this line. 
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Tabic 4: The isolated mean halo density profile for 
the mass < M >= 3 x 1O 12 /i _1 M . The symbols 
are the same as table 3. 



r/R vil 


< 8 > 


<5prob 


<<5'> 


c/ 

^prob 


1.0 


336 


323 


56.1 


52.4 


1.5 


119 


119 


15.1 


14.4 


2.0 


55.4 


65.6 


5.5 


4.2 


2.5 


30.1 


29.4 


3.1 


1.4 


3.0 


18.5 


18.2 


1.8 


0.36 


3.5 


12.3 


11.0 


1.4 


0.14 


4.0 


8.5 


7.0 


1.7 


-0.04 


4.5 


6.4 


6.0 


1.6 


0.03 


5.0 


5.2 


4.6 


2.5 


-0.14 


5.5 


4.5 


2.9 


1.8 


0.12 


6.0 


3.9 


2.7 


1.4 


-0.20 


6.5 


3.4 


2.0 


1.3 


-0.34 


7.0 


3.0 


1.5 


1.3 


-0.30 
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